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^Nj Abstract 

1—1 A permanent electric dipole moment (EDM) of a physical system requires time-reversal (T) and parity (P) violation. Experimental 

1 -^ programs are currently pushing the limits on EDMs in atoms, nuclei, and the neutron to regimes of fundamental theoretical 

I interest. Here we calculate the magnitude of the PT-violating EDM of ^He and the expected sensitivity of such a measurement 

fj to the underlying PT-violating interactions. Assuming that the coupling constants are of comparable magnitude for tt-, p-, and 

'—{ liJ-exchanges, we find that the pion-exchange contribution dominates. Our results suggest that a measurement of the He EDM 

^ IS complementary to the planned neutron and deuteron experiments, and could provide a powerful constraint for the theoretical 

'"^models of the pion-nucleon PT-violating interaction. 



^* 1. Introduction 
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' A permanent electric dipole moment (EDM) of a physical 
system would indicate direct violation of time-reversal (T) 
and parity (P) and thus CP violation through the CPT 
invariance. Presently there are several experimental pro- 
grams pushing the limits on EDMs in atoms, nuclei, and the 
neutron to regimes of fundamental interest. The Standard 
Model (SM) predicts values for the EDMs of these systems 
that are too small to be detected in the foreseeable future, 
and hence a measured nonzero EDM in any of these systems 
IS an unambiguous signal for a new source of CP violation 
and for physics beyond the SM. A new experimental scheme 
|1|2|3|4|5] for measuring EDMs of nuclei (stripped of their 
atomic electrons) in a magnetic storage ring suggests that 
the EDM of the deuteron could be measured to an accu- 
racy of better that 10~^^ e cm [4]. Unlike searches for CP- 
violating moments of the nucleus through measurements of 
atomic EDMs, a measurement for a stripped nucleus would 
not suffer from a suppression of the signal through atomic 
Schiff screening [6]. For this reason, the latter could repre- 
sent about an order of magnitude better sensitivity to the 
underlying CP- violating interaction than the present limit 
on the neutron EDM, dn [2]. Measurements using stripped 
nuclei in a magnetic storage ring are best suited to nuclei 
with small magnetic anomaly, making ■^He an ideal can- 



didate for a high precision measurement. Here we exam- 
ine the nuclear structure issues determining the EDM of 
•^He and calculate the matrix elements of the relevant op- 
erators using the no-core shell model [7] and Podolsky's 
method for implementing second-order perturbation the- 
ory [8]. An approximate and incomplete calculation for the 
■^He dipole exists in the literature [9], but here we present 
much more reliable calculations based on an exact solution 
of the three-body problem using several potential models 
for the nucleon-nucleon (NN) interaction, complemented 
with three-body forces. 

2. Sources of Nuclear i-',r violation 

A nuclear EDM consists of contributions from the fol- 
lowing sources: (i) the intrinsic EDMs of the proton and 
neutron, dp and d„; (ii) the polarization effect caused by 
the P-,r- violating {^"f) nuclear interaction, Hf,r^; (iii) the 
two-body -1^7^ meson-exchange charge operator appropriate 
for H^rp. 

The contribution due to nucleon EDMs, D^^^, which is 
purely one-body, can be easily evaluated by taking the ma- 
trix element 
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where {ip) is the nuclear state that has the maximal mag- 
netic quantum number. In the particular case of interest in 
this paper, \ip) = |0) is the ground state of ^He obtained 
by the diagonaHzation of the P,r-conserving interaction. 

In perturbation theory Hf,f induces a parity admixture 
to the nuclear state 
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where \n) are eigenstates of energy En and opposite parity 
from |0), which are calculated with the P-,r-conserving 
Hamiltonian. Hence, the polarization contribution £)^p°') 
can be simply calculated as 
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is the usual dipole operator projected in the z-direction. 

The contribution due to exchange charge, D'-'^^\ is typi- 
cally at the order of (w/c)^, and explicitly evaluated to be 
just a few percent of the polarization contribution for the 
deuteron case [lO]; we therefore ignore it and approximate 
the full two-body contribution, D^'^\ solely by the polar- 
ization term 
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Our calculation of the EDM of '^He therefore requires 
knowledge of both the individual EDMs of the nucleons 
and the f"t nuclear force. These very different quantities 
can only be related if some understanding exists of both 
the origin of the symmetry violation and its expression 
in strong-interaction observables. Constructing an effective 
field theory (EFT) that incorporates the symmetry viola- 
tion, as well as the dynamics underlying the usual strong- 
interaction physics in nucleons and nuclei, provides a suit- 
able framework. Chiral Perturbation Theory (xPT) supple- 
mented with a knowledge of the symmetry violation would 
be the appropriate EFT. To date only a single such calcula- 
tion exists [TT], and it was appHed to the one-nucleon sector, 
although further effort is underway |12I13) . The symmetry 
violation in that calculation was taken from the QCD 9 
term, which leads to an isoscalar f'T pion-nucleon interac- 
tion in leading order, unlike the most general case that in- 
cludes an isovector and an isotensor term, as well [H]. The 
non-analytic parts of the pion-loop diagrams |10jll|15fT6] 
that generate nucleon EDMs then provide an appropriate 
estimate of the EDMs of individual nucleons. These contri- 
butions are expected to dominate in the chiral limit |15j. 

In the absence of a xPT calculation of Hfr^ we revert to a 
conventional formulation in terms of a one- meson-exchange 
model. Including tt-, p-, and u;-meson exchanges, UJ the 
interaction is given by (see Refs. |10|17|18|19|20] ): 



^ Other mechanisms not included here are r;-exchange and two-7r- 
exchange, etc. 
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where to„ is the nucleon mass, G^ is defined as the product 
of a f't x-meson-nucleon coupHng gj (with T referring to 
the isospin) and its associated strong one, gxNN (e.g., G° = 
g-KNN 5° , where the interaction Lagrangian corresponding to 
these coupling constants is £ = N[ig.^,^til5 + 3°]t • TriV), 
yx{r) — e~™=" ''/(4 TT r) is the Yukawa function with a range 
determined by the mass of the exchanged x-meson, r = 
ri— r2, cr± = (Ti±cr2, and similarly for t± . Unless the sym- 
metry associated with the specific way that P, T violation 
is generated suppresses some of the couplings, one expects 
(by naturalness) that these fT meson-nucleon coupHngs 
are of similar magnitude, and this is roughly confirmed by a 
QCD sum rule calculation [2l].[2j We note, however, that in 
the (purely isoscalar) 9-term model of Ref . [11] the coupling 
constants G]^ and GJ vanish, and the coupling constants 
for the short-range operators are very small compared to 
the pion one [T3] . 

Because Hf,f violates parity and ^He is (largely) an S- 
wave nucleus, the matrix elements that define the EDM 
(see below) mostly involve 5"- to P-wave transitions. This 
has the combined effect of suppressing the short-range con- 
tributions and enhancing the long-range (i.e., pion) con- 
tribution, irrespective of the detailed nature of the force. 
Combined with the consideration that the short-range pa- 
rameters {Gp,uj) are not much larger than the pion ones, 
one can roughly expect the dominance of pion exchange. 

3. ^He in the ab initio No-Core Shell Model 

We solve the three-body problem in an ab initio no-core 
shell model (NCSM) framework ^. The ground-state wave 
function is obtained by a direct diagonalization of an effec- 
tive Hamiltonian in a truncated harmonic oscillator (HO) 
basis constructed in relative coordinates, as described in 
Ref. [22] ■ High-precision NN interactions, such as the local 
Argonne wis [23] and the non-local charge-dependent (CD) 
Bonn potential [24] interactions, are used to derive an ef- 
fective interaction in each model space via a unitary trans- 
formation [25] in a two-body cluster approximation. The 
Coulomb interaction between protons is also taken into ac- 
count. 

In addition to the phenomenological NN interaction 
models cited above, we consider two- and three-body inter- 
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actions derived from EFT. In a recent work [26j using the 
NCSM, the presently available NN potential at N^LO [22] 
and the three-nucleon (NNN) interaction at N^LO |28|29| 
have been applied to the calculation of various properties of 
s- and p-shell nuclei. In that study a preferred choice of the 
two NNN low-energy constants, cjj and ce, was found (and 
the fundamental importance of the chiral EFT NNN inter- 
action was demonstrated) by reproducing the structure of 
mid-p-shell nuclei. (Note that these interactions are fitted 
only for a momentum cutoff of 500 MeV, and therefore we 
are not able at this time to demonstrate a running of the 
observables with the cutoff.) This Hamiltonian was then 
used to calculate microscopically the photo-absorption 
cross section of ^He [30], while the full technical details on 
the local chiral EFT NNN interaction that was used were 
given in Ref. |3l]. We use an identical Hamiltonian in the 
present work, and we compare its predictions against the 
phenomenological potentials. 

In the NCSM the basis states are constructed using HO 
wave functions. Hence, all the calculations involve two pa- 
rameters: the HO frequency Q and N^ax, the number of os- 
cillator quanta included in the calculation. At large enough 
Nmax , the results become independent of the frequency, al- 
though the rate of convergence depends on fl. Thus, for 
short-range operators, one can expect a faster convergence 
for larger values of r2, as the characteristic length of the HO 
is 6 = 1/y/m^ Q,. The convergence also depends upon the 
P-,T-conserving interaction, Hq, used to solve the three- 
body problem. Thus the results obtained with Argonne wig 
show the slowest convergence, because the NN interaction 
has a more strongly repulsive core than the interactions 
obtained from EFT, which have faster convergence rates. 

The nucleonic contribution D^^> in Eq. JT]) involves only 
iJo, and is easily calculated once the three-body problem 
is solved. We therefore concentrate on the part involving 
Hff. 

Equation ((2]) suggests that in order to calculate the dipole 
moment one needs to obtain high-accuracy excited states 
of '^He in the continuum, an extremely difficult task in a 
NCSM framework, where the basis states are constructed 
using bound-state wave functions. The most straightfor- 
ward technique for evaluating Eq. ^ is to use Podolsky's 
method |8], in which |0) is obtained as the solution of the 
Schrodinger equation with an inhomogeneous term 

[E, -Hom^ Hfj |0) . (7) 

The exceptionally nice feature of this method is that con- 
tinuum states do not have to be explicitly calculated (they 
are, of course, implicitly included). In this sense the tech- 
nique is a relatively simple extension of bound-state meth- 
ods, which have been well studied and are robust. Moreover, 
in this approach the convergence of the EDM reduces to a 
large degree to the issue of the convergence of the ground 
state. 

We express the solution of Eq. (O as a superposition of 
a handful of vectors generated using the Lanczos algorithm 
[32] . Indeed, one can show that if we start with the inho- 



mogeneous part of Eq. ([7]) as the starting Lanczos vector 
|fi) = Hfr^ |0), the solution becomes [33] 
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where the summation over the index k runs over a finite and 
usually small number of iterations. The coefficients gk{E) 
are easily obtained using finite continued fractions [3l] . 

We alter this approach in practice for efficiency reasons. 
Because Eq. ([3]) is symmetrical in D^ and Hff, we are free 
to choose \vi) — Dz |0) as the starting vector. This allows 
us to isolate the two isospin contributions for Hff in each 
run. Once we compute a second vector, \v) = H*prp |0), the 
polarization contribution to the EDM is finally evaluated 
as 

^(Po0=2^gfc(i?o)(«k>. (9) 

k 

(We have verified that the altered approach gives the same 
results as the original one.) As a particular test case we 
have considered the electric polarizability 
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(where a is the fine structure constant), which reduces Eq. 
^ to aE = —2 a gi{Eo) (vi|ui). We estimate that the elec- 
tric polarizability of the "^He nucleus is 0.183 fm'^ for the 
Argonne vis potential, compared with 0.159 fm'^ reported 
in Ref. [35] for the same interaction. The 15% discrepancy 
is most likely the result of a difference in the theoretical 
approaches, as the result reported in Ref. [35] involves a 
matching of the ground-state energy to experiment (i.e., 
7.72 MeV), although the calculation gives 6.88 MeV bind- 
ing [3^ in the absence of three-body forces (our converged 
binding energy for Argonne vis is 6.92 MeV). Since the 
electric polarizability scales roughly with the inverse of the 
square of the binding energy, the discrepancy between the 
two results is reasonable. Moreover, we have made the ad- 
ditional check of the Levinger-Bethe sum rule [37], which in 
the case of tritium relates the total dipole strength to the 
charge radius, and we found it to be satisfied in all model 
spaces to a precision better than 10~^. Finally, the '^He po- 
larizability calculated using the two- and three-body chiral 
interactions is 0.148 fm'^, compared with 0.145 fm'^ with Ar- 
gonne Vis and Urbana IX [3^ two- and three-body forces. 
In both cases excellent agreement with the experimental 
binding energy is achieved. 

In a consistent approach the same transformation used 
to obtain the effective interaction should be used to con- 
struct the effective operators in truncated spaces. While 
this has been done in the past for general one- and two- 
body operators [38], such an approach is very cumbersome 
for the present investigation because both the dipole tran- 
sition operator and Hf,f change the parity of the states. 
We have therefore chosen not to renormalize the operators 
involved, except for the P-,T-conserving Hamiltonian. This 



problem is largely overcome by the fact that long-range op- 
erators (like the dipole) have been found to be insensitive 
to the renormalization in the two-body cluster approxima- 
tion [38j, which is the level of truncation for the effective 
interaction. We also point out that since Hff has short 
range, one can expect that the renormaHzation of this oper- 
ator would improve the convergence pattern, especially for 
small HO frequencies. As with all operators, the effect of 
the renormalization decreases as the size of the model space 
increases, so that in large model spaces (like the ones in the 
present calculation) this effect can be safely neglected and 
good convergence of observables is achieved. 

4. Results and Discussions 

We start the discussion of our results with the one-body 
contribution to the EDM of '^He. In Table [U we summarize 
the isoscalar and isovector contribution to Z?'^', which are 
decomposed into contributions proportional to their respec- 
tive coupHng constants {dp-\-dn for isoscalar, and dp — dn for 
isovector). All interactions give similar results, with only 
the Argonne wis result deviating more significantly from the 
others, albeit by less than 6%. We note that the coefficients 
in the upper and lower rows in Table[T]would be either 1/2 
or -1/2 if the nuclear forces between each pair of nucleons 
were taken to be equal (viz., the S'C/(4) limit, which implies 
that the neutron carries all of the nuclear spin) . 

Table 1 

The nucleonic contribution (in e fm) to the ^He EDM for different 
potential models. We decompose our results into contributions pro- 
portional to the nucleon isoscalar (dp + d„) and isovector {dp — dn) 
EDMs. 



converges to the experimental value already at N^ax ~ 22 
for most frequencies presented in the figure. 
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In Fig. [T] we present for four HO frequencies the running 
with N„iax of the EDM induced by the pion-exchange part 
of Hf,f. Two- and three-body EFT interactions have been 
used for this calculation, in order to obtain an accurate 
description of the ground-state of the three-body system. 
For the nuclear EDM we mix two types of operators: Hf,f, 
which is short range, and Dz, which is long range. The con- 
vergence pattern is therefore not as straightforward as pre- 
sented earlier in the discussion of convergence properties of 
general operators. The short-range part dominates the con- 
vergence pattern, and we thus observe faster convergence 
for larger frequencies (smaller HO parameter length) . This 
behavior is opposite to the convergence in the case of the 
electric polarizability, where we found faster convergence 
for smaller frequencies as expected for a long range opera- 
tor. Nevertheless, just as in Fig.[Tl the results become inde- 
pendent of the frequency at large Nmax- Note in the insert 
the convergence behavior of the ground-state energy, which 
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Figure 1. Isoscalar, isovector, and isotensor pion-exchange contribu- 
tions to the H fr^-'in&acedi EDM in ^He. We show four different fre- 
quencies in each case: f2 = 10 MeV (circles), f? = 20 MeV (squares), 
f! = 30 MeV (diamonds), and f! = 40 MeV (triangles). In the insert, 
we present the convergence of the ground-state energy, which in the 
limit of large Nmax approaches the experimental value (dashed line). 
Both NN and NNN EFT interactions have been used for diagonal- 
ization. 

Similar convergence patterns can be observed for the 
other meson exchanges as well as other potential models. 
In Table [2] we summarize these results. 

For pion exchange all potential models give basically the 
same result, as the long-range part (r > 1/m.^) of the 
■^He wave function shows negligible model dependence. It 
is interesting to note the effect of the three-body force by 
comparing the results with and without NNN interactions. 
When only the NN EFT interaction is used, the binding en- 
ergy is underestimated by about 500 keV. Therefore, since 
the ground-state energy is in the denominator of Eq. ([2]), 
one could naively expect that introducing the three-body 
forces (which increase the binding) decreases £)(p°'). In- 



Table 2 

The polarization contribution to '^He EDM (in units of e fm), decomposed as coefficients of ff^ , where x stands for tt, p, or lo meson exchanges. 
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stead we obtain nearly the same results for both isoscalar 
and isovector contributions. This implies that the NNN 
interaction reshuffles the strength to compensate for the 
change in binding energy, most likely at low energies. This 
is not a surprise, because it was already found previously 
that the main effect of the three-body forces for the dipole 
response is an attenuation of the peak region at low ener- 
gies both in the three- [35] and four-body [30] systems. 

In the isoscalar pion-exchange channel our result (0.015) 
is about 50% larger than an existing work [9] , which yielded 
0.010. Since the Reid soft core NN potential used in Ref. 
[9] has a much more repulsive core than even Argonne vi^ , 
that isoscalar contribution to the EDM is in line with our 
findings. Moreover, the approximate solution of the three- 
body problem (as opposed to the current work, in which the 
calculation is exact) can induce uncontrolled uncertainties. 
Finally, their calculation of the isovector term did not ex- 
haust all the possible spin-isospin combinations, while the 
isotensor contribution was not computed. 

For p- and w-exchanges one immediately sees that their 
corresponding coefficients are at most 10% of the pion- 
exchange ones, because only the short-range wave function 
(r < l/nip^^) contributes substantially. Sensitivity of the 
short-range f'f potentials to the short-range model de- 
pendence of the wave functions produces matrix-element 
variations as large as 50% for some channels. While a de- 
tailed explanation for the model dependence is too intri- 
cate to be disentangled, one can roughly see the general 
trend that the calculation using Argonne wis gives consis- 
tently smaller results than ones using CD Bonn and chiral 
EFT, as Argonne wig has a harder core than the other two. 
The behavior of short-range f't nuclear potentials in xPT 
strongly suggests that the short-range coupling constants 
are no larger than the pion ones, and may be significantly 
smaller. If this holds, pion-exchange will produce the domi- 
nant contribution to nuclear EDMs, with the heavy-meson 
exchanges (the short-range interaction) giving roughly a 
10% correction (or less) to the pion-exchange contribution. 
This suppression due to P-wave intermediate nuclear states 
was discussed earlier and is in accord with calculations in 
heavier systems [20j . 

Assuming the dominance of pion exchange, Z?'^' has an 
almost model-independent expression 

i:»(2) r^ (0.015 G^ + 0.023 G^ + 0.036 G2)efm. (11) 

The single-nucleon EDMs can be estimated using the non- 



analytic term that results from calculating the one-pion 
loop diagram, which dominates in the chiral Hmit (see, for 
example, Refs. [10|ll|14|15|l6] l 
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where e is the proton charge. Folding this result into 
D'-^'> and using the physical nucleon and pion masses 
(ln(77T,jv/m7r) ~ 1.90) we get 

ZJfi' w 0.009 (G°-G^)efm . (13) 

The total EDM of '^He is therefore estimated to be 

= (0.024 G° + 0.023 G^ + 0.027(5^) e fm . (14) 

Calculating the EDMs of the neutron and deuteron [lO] 
using Eq. I|12p and assuming pion-exchange dominance, one 
can see from Table [3] that an EDM measurement in '^He 
is complementary to the former two. Assuming that sim- 
ilar sensitivities can be reached in these three measure- 
ments, the I^T' pion- nucleon coupHng constants could be 
well-constrained if the assumption of pion-exchange domi- 
nance holds. 

Table 3 

Complementarity of the ^He EDM to the neutron and deuteron 

EDMs. We present the theoretical estimations (in units of e fm) 

of neutron, deuteron, and He decomposed into contributions pro- 

— 12 

portional to GV ' , while assuming the dominance of pion-exchange 
forces in H^r^ and estimating nucleon EDMs via pion loops. 
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5. Summary 

In this paper, we have calculated the nuclear EDM of ''He, 
which arises from the intrinsic EDMs of nucleons and the P- 
,T- violating nucleon-nucleon interaction. Several potential 
models for the P-,T-conserving nuclear interaction (includ- 
ing the latest-generation NN and NNN chiral EFT forces) 
have been used in order to obtain the solution to the nuclear 
three-body problem. The results obtained with these po- 
tential models agree within 25% in the |^7^ pion-exchange 



sector. Though larger spreads in f't p- and w-exchanges are 
found (as the results sensitively depend on the wave func- 
tions at short range), we expect them to be non-essential 
as pion-exchange completely dominates the observable (un- 
less the f't parameters associated with heavy-meson ex- 
changes are significantly larger than the ones for pion ex- 
change, which is not expected). We further demonstrate 
that a measurement of the '^He EDM would be complemen- 
tary to those of the neutron and deuteron, and in combi- 
nation they can be used to put stringent constraints on the 
three P-,r- violating pion-nucleon coupling constants. We 
therefore strongly encourage experimentalists to consider 
such a '^He measurement in a storage ring, in addition to 
the existing deuteron proposal [4]. 
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